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ABSTRACT 

We examine the age- dependence of dark matter halo clustering in an unprece- 
dented accuracy using a set of 7 high-resolution cosmological simulations each 
with N = 1024 3 particles. We measure the bias parameters for halos over a large 
mass range using the cross-power-spectrum method that can effectively suppress 
the random noise even in the sparse sampling of the most massive halos. This 



enables us to find, for the first time, that younger halos are more strongly clus- 
tered than older ones for halo masses M > 10M*, where M* is the characteristic 
nonlinear mass scale. For M < M*, our results confirm the previous finding of 
Gao et al. that older halos are clustered more strongly than the younger ones. 
We also study the halo bias as a function of halo concentration, and find that 
the concentration dependence is weaker than the age dependence for M < M*, 
but stronger for M > 50M*. The accurate and robust measurement of the age 
dependences of halo bias points to a limitation of the simple excursion set the- 
ory which predicts that the formation and structure of a halo of given mass is 
independent of its environment. 



Subject headings: gravitational lensing — dark matter — cosmology: theory — 
galaxies: formation 
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Introduction 



The formation of dark matter halos plays a central role in the studies of galaxy for- 
mation as well as of the large-scale structure of the universe. A widely used analytical 
theory for the formation of hal os is the extended Press-Schechter (PS thereafter) formalism 



(IBond et al 



1991 



Bower! 1 199 if) , which can be used to model the mass functio n of dark ha- 



los ([Press fc Schechterl 11974 ). the bias parameter as a function of halo mass (IMo fc White 



1996, hereafter MW96), as well as the formati on histories o f dark halos (ILacey fc Cole 



1993 



1994 ). In the simplest excursion-set model of IBond et al.l (1199 ll ). the physical properties of 
dark halos are expected to depend only on halo mass, but not on large-scale environments. 
However, it has been known for almost a decade that neither the PS mass function nor 
the MW96 bias model matches well iV-body results for halo masses M < M*, where M* is 
defined such that the rms linear density fluctuation within a sphere of mass i s S n = 1.68 



Ming 



1998 



Lee fc Shandarin 



1998 



Sheth fc Lemsonl Il999l : IPorciani et al.lll999l : Ijindll999 



Sheth fc Tormenj|l999l ). ISheth et al.l (120011 ) proposed an ellipsoid collapse model to replace 
the spherical collapse model in the PS theory, and found a better agreement between the 
theory and TV-body simulations in both the mass function and the bias function. Earlier 
investigations with cosmological iV-body simulations did not find any strong environmental 
dependence of halo propertie s , such as the spin pa rameter, concentration and formation time 
( Lemson fc Kauffmannlll999l ; IPercival et al.l 120031 ). However, the poor resolutions and small 
dynamical ranges covered in these simulations make it difficult to detect any signals that are 
relatively weak or outside the dynamical range, and so a significant age-dependence of halo 
clustering cannot be ruled out. Note also that there we re earlier attempts to develop empir- 



ical models for the age-depe ndence of halo clustering (ITaruya &: Sutol |2000| ; lHamana et al. 
200ll ; lYoshikawa et al.llioOlh . 



With the use of a very large iV-body simulation, iGao et al.l (120051 ) recently demonstrated 
convincingly that the clustering strength of dark halos does depend on halo formation time. 
In particular, they found that this dependence is strong o nly for halos w it h M -C M* and 
becomes very weak for M > M* . This also explains why IPercival et al.l (120031 ) could not 
detect an age dep endence, because they focused on halos with M > M* at high redshifts. 



Wang et al.l (120061 ) have examined the physical process that may be responsible for the age 
dependence of halo clustering. They found that halos embedded in dense environments 
accrete mass less efficiently than the spherical collapse model predicts, because the matter 
to be a ccreted is 'heated' by the large-scale structure (like the pancake heating considered 
in, e.g.. lMo et al.ll2005l ; lLedl2006l ). This explains why the old population of small halos has a 
higher bias than the young population of the same mass. It also qualitatively explains why 
the PS mass function and bias functions deviate from iV-body simulations at small masses. 
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It is well known that the concentration parameter of halos at a given mass, c(M), has a 
broad log-normal distribution, and is correlated with the halo for mation epoch Zf mm in such 



a way that halos of earlier formation have a higher concentrati on (jjing 2000 ; Bullock et al. 



200ll ; |jing fc Sutolbopj fwechsler et aPhooi Ehao et alil2003al lb}). Therefore the age depen- 
dence of halo bias should yield a concentration dep endence of halo bias. Such dependence 
is indeed observed in recent iV-body simulations. Wechsler et al] ( 2006 ) found that, for 
masses M < M*, halos with higher concentrations are more strongly clustered, as expected 
from the age dependence and the correla t ion b etween the age and concentration of dark 
halos. More interestingly, I Wechsler et al.l (120061 ) found reversed concentration dependence 
for M > M*, in the sense t hat halos with hi gher concentrations are actually less biased. 
This result was not seen by iGao et al.l ( 2005 ) who explored the age dependence only for 
M < 10M*. More recently, IWetzel et al.l (120061 ) examined both the age dependence and con- 
centration dependence of halo cluste ring for massive halos . While they found concentration 
dependence similar to that found by I Wechsler et al.l (120061). they did not detect any signifi- 
cant dependence on formation epoch. Note, however, that IWechsler et al.l (120061 ) used halos 
identified at different redshi fts to increase the dynam ical range probed by their simulations, 
while both lGao et al.l (120051 ) and lWetzel et al.l (120061 ) used halos identified at the same time. 
It is unclear how to make a detailed comparison between the different results. 

The environmental dependence of halo formation and structure has important implica- 
tions not only for improving the PS theory but also for improving semi-analytical models of 
galaxy formation based on PS merger trees and current halo occupation distribution (HOD) 
models. In this Letter, we use a large set of cosmological simulations, each with 1024 3 par- 
ticles, to investigate in detail the concentration and formation epoch dependence of halo 
clustering. 



Simulations 



The model considered here is a canonical spatially-flat Cold Dark Matter (CDM) model 
with the density parameter Q m = 0.268, the cosmological constant = 0.732, the Hubble 
constant h = 0.71, and the baryon density parameter Qf, = 0.045. The primordial density 
field is assumed to be Gaussian with a scale-invar iant power spectrum oc k. For the linear 
spectrum, we adopt the fitting transfer function of lEisenstein fc Hul (119981 ). and the normal- 
ization is set by erg = 0.85, where <7g is the present linear RMS density fluctuation within a 
sphere of radius 8 h~ 1 Mpc. 



We use an upgraded version of the Particle-Particle-Particle- Mesh (P 3 M) code of lJing fc Suto 



( 119981 . |2002| ) to simulate structure formation in the universe. The code has now incorporated 
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the multiple level P 3 M gravity solver for high density regions (jJing &: Sutol 120001 ) . In order 
to have a large mass resolution range, we run a total of 7 simulations with 1024 3 particles in 
different simulation boxes (Tabled]). The simulations were run on an SGI Altix 350 with 16 
processors with OPENMP parallelization in Shanghai Astronomical Observatory. We have 
4 realizations with boxsize 1800 /i _1 Mpc in order to reliably measure the bias of the most 
massive halos, which is the focus of the present Letter. 

Dark matter halos are identified using the standard Friends-of-Friends algorithm with 
a linking length b equal to 0.2 times the mean particle separation. Unbound particles (with 
positive binding energy) are excluded. We use halos containing 100 particles or more, and 
our analysis covers halos with masses ranging from 2 x 10 11 /i _1 M Q to 10 15 h~ 1 M Q . 



3. Clustering of dark matter halos 



3.1. Cross power spectrum and the halo bias 

Once the number density field of dark matter halos, n/j(r), and the density field of dark 
matter p m (r), are given, the conventional way to estimate the halo bias factor b is to use the 
definition, b = (£/ih(r)/£ mm (r)) 1 ' /2 , where £hh{r) is the two-point correlation function of dark 
halos, and C, mm (r) that of da rk matter. H ere we adopt a slightly different approach by using 
the cross power spectrum of iJina (119991 ) . We first Fourier transform both ri^(r) and p m (r) 
into rih(k.) and p m (k), and then measure the bias factor through 

Phm(k) p m (n h (k)p m (k)) 



b = 



Pmm(k) U h (p m (k)p m (k)) 



1) 



As shown by I Jingj ( 119991 ). this method suppresses Poisson noise due to a limited number 
of dark halos more effectively than the correlation method, because the number of dark 
particles is in general much larger than that of dark halos. This is particularly important 
for the present study, as we need to accurately determine the bias for massive halos whose 
number density is low. We have estimated the bias factors for the four realizations of 
the 1800 /i _1 Mpc box simulations, using both the correlation method and the cross power 
spectrum method. The average values of the bias over the four realizations obtained with 
the two methods are consistent, but the scatter with the correlation method is about a factor 
of 2 larger than that with the cross power spectrum. The cross power spectrum method is 
therefore preferred. 



Errors on b(k) are estimated following I Jingi (119991 ) : ten random samples are generated 
for each simulated halo sample, and the scatter of b(k) among them is used as the error in 
b(k) for the real sample. We have compared the error so obtained with that estimated from 
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the scatter among the 4 realizations of the 1800 h 1 Mpc box simulation. The errors given 
by the two methods are comparable. 



3.2. Dependence on the halo formation epoch 

Our analysis focuses on the halo populations at the present time, z = 0. Halos are 
divided into mass bins (specified by [777,1,777.2] with 777 2 = 2mi) according to their masses. 
For each mass bin, we sort the halos into 5 populations according to formation redshift, Zf, 
which is defined as the redshift at which the mass of the most massive progenitor of a halo is 
equal to half of its mass. Although analysis was carried out for all the five populations, here 
we present the results only for the 20% youngest halos with the smallest Zf and for the 20% 
oldest ones with the largest Zf. In what follows, these two populations are simply referred to 
as the young and the old populations, respectively. To ensure that the bias obtained here is in 
the linear regime, we only consider Fourier modes with k < k max , where k max = 0.09 /iMpc -1 
is a wavenumber at which the variance A 2 (k max ) = kf nax P mm (k rnax )/2ir 2 = 0.25. Figured] 
plots examples of the bias factor for massive halos with M = 35M* and M = 134M*, where 
M* = 5 x 10 12 /i _1 M Q , in the 1800 /i _1 Mpc box simulations. The top-left panel shows the 
results for the young population, while the top-right panel shows the same measurement for 
the old population. The results clearly show that the bias factors are nearly scale independent 
in the linear regime. The solid line in each panel is the mean value obtained by a least square 
fit to the data points at k < k max . The bias factor for the young population is about 10% 
higher than that for old population. 

We present the mean bias factor b(M) as a function of mass obtained from different 
simulations in the top panel of Figure [2J For a given halo mass, the measurement of b(M) 
is more accurate in a simulation of a larger volume. Therefore we plot the data points for 
each simulation only up to a halo mass which corresponds to 100 particles in the next larger 
simulation, except for the 1800 /i _1 Mpc box simulation, where we plot all the data points 
available. Note that the bias factors obtained from different simulations agree extremely 
well with each other, and so the set of simulations used here allows us to explore the age 
dependence of halo bias over four orders of magnitude in halo mass. For M < 3M*, the 
old population (red symbols) is more strongly biased (has stronger clustering) than the 
young popu lation (blue symbols), which is in a good quantitative agreement with the results 
obtained bv lGao et al. hood, cyan lines). With the Millennium Simulation with a boxsize of 



500 ft, -1 Mpc, iGao et al.l (120051 ) was not able to explore accurately the age dependence of halo 
clustering for higher halo mass. Our results clearly show that the difference between the old 
and young populations decreases with increasing halo mass up to M ~ 10M*, at which the 
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young population starts to surpass the old population in the bias factor. For M > 20M*, the 
bias factor for the young population is about 10% higher than that of the old population, 
with weak dependence on halo mass. Although the difference in the bias factor between the 
old and young populations is small at the high mass end, it is detected at a high statistical 
confidence (~ 10<t). 



3.3. Dependence on halo concentration 



The concentration of each halo is obtained following the fitting method of iJing) (120001 ) . 
The density distribution within each halo is fitted with a NFW profile to obtain the scale 
radius r s , and the concentration is defined as c = r v /r s , where r v is the virial radius within 
which the mean density is 361 times the mean density of the universe. Here only halos with 
320 particles or more are used, because the concentr ation may not be me asured accurately 
for halos containing smaller number of particles (e.g. IWechsler et al.ll2006l ). 



The lower two panels of Figure [T] show the bias factor as a function of wavenumber for 
massive halos with M = 35M* and M = 134M*. Results for the least concentrated 20% 
are plotted in the lower left panel, while those for the most concentrated 20% are in the 
lower right panel. Here again, the bias factor is almost scale-independent. The amplitude 
of the bias factor for massive halos clearly depends on concentration, with halos with higher 
concentration less strongly biased. To see how the concentration-dependence changes with 
halo mass, we plot in the lower panel of Figure [2]the bias factor for the most concentraed 20% 
and the least concentrated 20% of the halos in each of the mass bins. Note that there is good 
agreement between simulations with different box sizes, suggesting that 320 particles may 
be sufficient to sample the concentration for the purpose of the present paper. Our results 
show clearly that the more concentrated halos have a larger bias for M < M*, but the tren d 
is reversed for M > M#, in qualitative agreement with the results in lWechsler et al.l (120061 ) . 



Comparing the results here with the dependence on the formation epoch (the upper 
panel), we see that the concentration dependence is weaker than the age dependence for 
M < M#, but stronger for M 3> M*. The difference in the bias factor between the most 
concentrated 20% and the least concentrated 20% is about 25% for M = 10M* - 100M*, larger 
than the 10% between the youngest 20% and oldest 20%. This may be why concentration 
dependence wa s but age dependence was not f o und for halos with M > in previou s 



investigations (IGao et al.l l2005t IZhu et al.l 12009 : IWechsler et al.l 120061 : IWetzel et all |2006|). 
Our results also show that the concentration dependence reverses almost exactly at M = M*, 
while the reversal of age dependence occurs at a much larger mass, M ~ 10M*. Finally, 
the concentration dependence does not seem to become stronger with increasing mass for 
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M > M*. In fact, the difference in the bias factor between the most concentrated 20% and 
the least concentrated 20% is only about 10% at 250M*. 



3.4. The effect of finite mass bins 

So far we have used a finite mass bin [Mx, M 2 ] with M 2 = 2M\ to study the age and 
concentration dependence of the bias for halos of a given mass. Since more massive halos 
on average have younger ages and smaller concentration, and since the bias factor increases 
strongly with halo mass at the high mass end, the use of a finite mass bin to represent a givn 
mass may artificially introduce age and concentration dependence. To check thie effect, we 
repeat our analysis using narrower mass bins with M 2 = 1.2M\. The results are plotted as 
the solid and dashed lines in Figure [2j The results are almost indistinguishable from those 
obtained with M 2 = 2Mi, demonstrating that the mass bins we used is sufficiently small. 



4. Summary and Discussion 

Our major findings in this paper can be summarized as follows: 

(i) The younger halo population is significantly less clustered than the old one for M < 

10M,; 

(ii) For M > 10M*, the age dependence is reversed, and the bias factor of the youngest 20% 

of the halos with a given mass is approximately 10% higher than that of the oldest 
20% of the same mass; 

(iii) When the halos are divided into subsamples according to concentration parameter, 
the halo bias is larger (smaller) for halos with higher concentrations for M < M* 
(M > M*) . 



The first result is a confirmation of the result of iGao et al.l (120051 ) with the use of a much 
larger data set and a different clustering measure (our cross power spectrum versus their 
auto-correlation function). Our use of a large set of simulations in large boxes enables us 
to quantify, for the first time, the age dependence of halo bias for M > 10M*. The weak 
age dependence for massive halos relative to low-mass halos ones may re flect the fact that 



these halos have a narrower distribution in formation time (see Fig.l of iKitayama &: Suto 



1996bl ). but the opposite trend, although weak , is no t easy to explain. The third result is 



qualitatively in agreement with IWechsler et al.l (120061 ). However, their results are based on 
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halos identified at different times, while ours are based on halos identified at the same tme 
as it should be. Our results also show that the concentration dependence is weaker than the 
age dependence for M < M* but stronger for M ~ 50M*. There is also indication that both 
the age and concentration dependence becomes weaker at the very massive end. 

The difference in the age and concentration dependence implie s that caution m ust be 
taken in comparing model predictions with observational results (e.g. lYang et al.ll2006l ). since 
it is unclear whether the formation epoch, as defined by Zf, or concentration, or even some 
other properties of a halo is more important in determining the properties of the galaxies 
that form in it. 

The dependence of halo bias on format ion epoch and concentra t ion implies that the sim- 



ple ex cursion set theory of halo formation ( Press fc Schechterlll974l ; lBowerlll99ll ; iBond et al. 
199fl ) is not accurate. As discussed in I Wang et al.l (120061 ). the problem with the simple 



excursion set theory is that spherical collapse model, which neglects large-scale tidal field, 
over-predicts the collapse of small halos (<C M*) in high density environments. They ar- 
gue, however, that this dynamical effect should be smaller for large halos (M ^> M#). It is 
possible that the large-scale tidal field also plays a role in the formation of massive halos. 
But instead of truncating mass accretion, as is for the case of low-mass halos, the large-scale 
tidal field may delay the accretion, enhancing the accretion at later times. This possibility 
should be checked further by examining in detail the accretion histories and environments 
of massive halos in simulations. 



The extended PS theory (ILacey fc Cold Il993l . Il994 ; ISasakil Il994 ; iKitayama fc Suto 



1996af) has been implemented in a variety of cosmolog i cal studies with clus ters of galax- 
ies ( iKitayama fc Sutolll996bl . 119971 ; iTaruya fc Sutoll2000l ; lHamana et al.ll200ll ). The present 
results for massive halos point to a limitation of the theory in such applications, and indeed 
should prompt the exploration of an improved theoretical framework in the future. 
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Table 1. Simulation parameters 



boxsize( h 1 Mpc) 


particles 


realizations 


^particle 


300 


1024 3 


1 


1.8 x 1O 9 /i _1 M 


600 


1024 3 


1 


1.5 x 10 10 h- 1 M Q 


1200 


1024 3 


1 


1.2 x 10 11 h- 1 M Q 


1800 


1024 3 


4 


4.0 x 1O 11 /i~ 1 M 
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Fig. 1. — The ratio of the cross power spectrum P/ im (fc) of halos with background dark 
matter to the matter power spectrum P mm (k) versus the wavelength in the linear regime. 
The top left and right panels are for the 20% youngest and 20% oldest halos respectively, 
where the age is defined as the formation redshift Zf. The bottom left and right panels show 
the results for the 20% halos that have the highest and lowest concentrations respectively. 
The mass of the halos, in unit of the characteristic mass, is given in the panels. 
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Fig. 2. — The bias factor of dark matter halos as a function of the halo mass. In the upper 
panel, we show the dependence on the halo formation epoch, with the symbols in red being 
for the 20% o ldest halos and th ose in blue being for 20% youngest halos. The cyan lines are 
the results of iGao et al.l (120051 ) for comparison. The lower panel shows the dependence on 
the halo concentration, with the red and blue colors being for those with the 20% highest 
and lowest concentrations respectively. The open triangles, open circles, filled triangles, and 
filled circles are from the simulations of boxsize 300 h~ l Mpc, 600 h~ 1 Mpc, 1200 h~ 1 Mpc and 
1800 /i _1 Mpc respectively. The solid and dashed lines are not the lines connecting the data 
points, but are for the results of young (or less concentrated) and old (or more concentrated) 
halos respectively estimated for a very narrow mass bin M ± 0.1M. 



